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In this study, a two-state mechanochemical model is presented to describe the dynamic insta- 
bihty of microtubules (MTs) in cells. The MTs switches between two states, assembly state and 
disassembly state. In assembly state, the growth of MTs includes two processes: free GTP-tubulin 
binding to the tip of protofilament (PF) and conformation change of PF, during which the first 
tubulin unit which curls outwards is rearranged into MT surface using the energy released from the 
hydrolysis of GTP in the penultimate tubulin unit. In disassembly state, the shortening of MTs 
includes also two processes, the release of GDP-tibulin from the tip of PF and one new tubulin unit 
curls out of the MT surface. Switches between these two states, which are usually called rescue and 
catastrophe, happen stochastically with external force dependent rates. Using this two-state model 
with parameters obtained by fitting the recent experimental data, detailed properties of MT growth 
are obtained, we find that MT is mainly in assembly state, its mean growth velocity increases with 
external force and GTP-tubulin concentration, MT will shorten in average without external force. 
To know more about the external force and GTP-tubulin concentration dependent properties of 
MT growth, and for the sake of the future experimental verification of this two-state model, eleven 
critical forces are defined and numerically discussed. 



I. INTRODUCTION 

In eukaryotic cells, microtubules (MTs) serve as tracks 
for motor proteins [l|-|6|, give shape to cells, and form 
rigid cores of organelles [THlOj. They also pla y essen- 



tial roles in the chromosome segregation |lll - ll8l | . During 
cell division, MTs in spindle constantly grow and shorten 
by addition and loss of enzyme tubulin (GTPase) from 
their tips, then the attached duplicated chromosomes are 
stretched apart (through two kinetochores) from one an- 
other by the opposing forces (produced by MTs based 
on different spindles). Recently, many theoretical mod- 
els have been designed to understand the roles of MTs 
during chromosome segregation |19l425| . One essential 
point to understand how MTs help chromosome segrega- 
tion during cell division is to know the mechanism of MT 
growth and shortening. In this study, inspire d by the 
mechanochemical model for molecular motors 12611 . the 



GTP-cap model and catch bonds model for MT |8l.l27l|. 
a two-state mechanochemical model will be presented. 

Electron microscopy indicates MT is composed of n 
parallel protofilaments (PFs, usually 12 < n < 15 and 
n = 13 is used in this study) which form a hollow cylinder 
0,BEil- Each PF is a filament that made of head-to-tail 



associated a/S heterodimers. At the tip of MT, PFs curl 
outward from the MT cylinder surface. The tip might be 
in shrinking GDP-cap state or growing GTP-cap state. 
In contrast to the tip in shrinking GDP state, the growing 
GTP tip is fairly straight. Or in other words, in GTP-cap 
state, the angle between the curled out segment of PFs 
and MT surface is less than that in GDP state. In this 
study, we will only consider the growth and shortening 
of one single PF, and assume that each step of growth 
and shortening of one PF contributes to L (nm) of the 
growth and shortening of the whole MT. Intuitively, L = 
Li/n with Li the length of one af3 hetrodimer. In the 
numerical calculations, L = 8 nm/13 «0.615 nm is used 
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Our two-state mechanochemical model for PF growth 
and shortening is schematically depicted in Fig. [TJa), and 
mathematically described by a two-line Markov chain in 
Fig. [IJb). In this model, PF stochastically switches be- 
tween two states: assembly state and disassembly state. 
During assembly state, PF grows through two processes, 
(i) 1-^2: free GTP-tubuhn binding process with GTP- 
tubulin concentration (denoted by [Tubulin]) dependent 
rate constant ki, and (ii) 2^1: PF conformation 
change process, during which, using energy released from 
GTP hydrolysis, the curled PF segment is straightened 
with one PF unit (i.e. one a/3 heterodimer) rearranged 
into the MT surface, i.e. to parallel the MT axis ap- 
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proximately. During disassembly state, each step of PF 
shortening includes also two processes, (i') 2' ^ 1': dis- 
association of GDP-tubulin from PF tip to environment 
and (ii') 1' 2': one new PF unit curls out from the 
MX surface (during which phosphate is released from the 
tip tubulin unit simultaneously). 

The two-state model presented here can be regarded 
as a generalization of the one employed by B. Akiyoshi 
et al to explain their experimental data [27], which can 
be depicted by Fig. [2ja) [our corresponding generalized 
two-state model including bead detachment from MT is 
depicted in Fig. [2l[b), see Sec. II. A for detailed discus- 
sion]. The reasons that we prefer to use this general- 
ized model are that, the simple model of B. Akiyoshi et 
al cannot fit the measured attachment lifetime of bead 
on MT well (see Fig. 4(a) in 0j), and moreover, the 
measurements in 301-1321] indicate that the rate of catas- 
trophe, i.e. transition from elongation to shortening, de- 
pendent on GTP-tubulin concentration of the solution. 
However, for the simple model depicted in Fig. [2Ia), the 
catastrophe rate kc is independent of GTP-tubulin con- 
centration (it is biochemically reasonable to assume that 
the elongation rate fci depend on GTP-tubulin concentra- 
tion, fci = k\ [Tubulin] , but with no reasons to write kc as 
as a function of [Tubulin]). We will see from the Results 
section that, for our generalized model, the catastrophe 
rate does change with [Tubulin] , since GTP-tubulin con- 
centration will change the probabilities of PF in states 1 
and 2, and consequently change the transition rate from 
assembly state to disassembly state. At the same time, 
for the simple model depicted in Fig. [2Ia), the distribu- 
tion of catastrophe time is an exponential. However, the 
experimental measurement under a particular situation 
indicates this distribution is clearly not an exponential 
[3^ 44]. It should be pointed out, although our model 
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presented here looks more complex, there are only two 
more parameters than the one depicted in Fig. [D^a) [4^. 



The organization of this paper is as follows. The two- 
state mechanochemical model will be presented and theo- 
retically studied in the next section, and then in Sec. HI, 
based on the model parameters obtained by fitting the 
experimental data mainly obtained in f27|, properties of 
MT growth and shortening are numerically studied, in- 
cluding its external force and GTP-tubulin concentration 
dependent growth and shortening speeds, mean dwell 
times in assembly and disassembly state, mean growth or 
shortening length before the bead, used in experiments 
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FIG. 1: Schematic depiction of the two-state mechanochem- 
ical model of protofilament (PF) growth and shortening (a) 
and its corresponding two-line Markov chain (b). In assem- 
bly state, the growth of PF accomplished by two processes, 
one GTP-tubulin binds to the tip of PF (with GTP-tubulin 
concentration dependent rate ki) and one PF unit rearranges 
into the MT surface (with external force dependent rate ^2). 
The energy used in the second process comes from the GTP 
hydrolysis in the penultimate tubulin unit. One tubulin unit 
binding to the tip of PF is assumed to be equivalent to L 
(nm) growth of the whole MT (L = 0.615 nm is used in this 
study [1^, [2^). Similarly, in disassembly state, the shortening 
of PF also includes two processes, one PF unit detaches from 
the tip of PF and one new PF unit curls out the MT surface. 
In this depiction, the same as in [2^, a segment 5 dimers in 
length is assumed to curls out from the MT surface. 



to apply external force, detachment from MT. To know 
more properties about the MT dynamics, eleven critical 
forces (detailed definitions will be given in Sec. HI) are 
also numerically discussed in Sec. HI. Finally, Sec IV 
includes conclusions and remarks. 



II. TWO-STATE MECHANOCHEMICAL 
MODEL OF PROTOFILAMENT 



As the schematic depiction in Fig. [U PF might be 
in two states, assembly (growth) state and disassem- 
bly (shortening) state. Each of the two states includes 
two sub-states, denoted by 1, 2 and 1', 2' respectively. 
Let pi,P2 be the probabilities that PF in assembly sub- 
states 1 and 2 respectively, and pi , p2 be the probabili- 
ties that the PF in disassembly sub-states 1' and 2', then 
pi,P2,Pi,P2 are governed by the following master equa- 
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tion 



dpi/dt =k2P2 - kiPi + Kpi - kcPi, 
dp2/dt =fcipi - /C2P2, 
dpi/dt ^kiP2 - kspi - krpi + kcPi, 
dp2/dt =k-ipi - kiP2. 



(1) 



Where fci is the rate of GTP-tubuhn bmding to the tip 
of PF, ^2 is the rate of PF reahgnment with one new 
unit lying into the MT surface, k^ is the dissociation rate 
of GDP-tubuhn from the tip of PF, and ki is the rate of 
curUng out of one tubuhn unit from the MT surface (with 
Pi release simultaneously). The steady state solution of 
Eq. © is 



Pi 



P2 



k2 kr 



k4 



ki 



k?. kr 



(2) 
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One can easily show that the mean steady state velocity 
of MT growth or shortening is [33|, |34 1 



V = {k2P2 - kip2)L = (fci - k3kc/kr)piL, (3) 

where L is the effective step size of MT growth corre- 
sponding to one step growth of one PF, and V < means 
MT is shortening in long time average with speed ^V. 

Let pi,P2 be the probabilities that PF in sub-state 1 
and sub-state 2 respectively, provided the PF is in as- 
sembly state, then pi,p2 satisfy 



dpi/dt =k2P2 - kipi = -dp2/dt. 



(4) 



One can easily get that, at steady state, the mean growth 
speed of MT with a PF in assembly state is 



Vg = k2P2L 



kik2L 
ki + k2' 



(5) 



Similarly, the mean shortening speed of MT with a PF 
in disassembly state is 



k^k^L 
fca -I- k4 



(6) 



Modified model according to experiments: 
including bead detachment from MT 



To know the model parameters ki 



1, 



, 4 and 



kcykr-i we need to fit the model with experimental data. 
In recent experiments Akiyoshi et al attached a bead 
prepared with kinetochore particles to the growing end of 
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FIG. 2: (a) Schematic depiction of the two-state model used 
by B. Akiyoshi et al in [271 ]. In which, both the assembly 
and disassembly of PF are assumed to include only one pro- 
cess, described by rates fci and fe2 respectively, (b) Modified 
mechanochemical model with bead detachment. In the ex- 
periments of _27], Akiyoshi et al attached a bead prepared 
with kinetochore to the growing end of MTs, and measured 
not only the force dependent mean growing and shortening 
speeds, switch rates between assembly and disassembly states 
(i.e. rates of rescue and catastrophe), but also the mean life- 
time of the bead on MTs, and the rates of bead detachment 
during assembly and disassembly states respectively. There- 
fore, to get the model parameters and know more properties 
of MT growth and shortening, this modified model is used 
in this study . The main difference between these two mod- 
els is that, the rate of detachment from assembly state and 
the rate of catastrophe in model (b) depend on GTP-tubulin 
concentration [Tubulin], but they do not in model (a). 



MTs, and constant tension was applied to bead using a 
servo-controlled laser trap. In their experiments, not only 
the force dependent mean growth and shortening speeds 
of MTs, the rates of rescue and catastrophe, but also the 
force dependent mean lifetime, during which the bead is 
keeping attachment to MT, and mean detachment rates 
of the bead during assembly and disassembly states are 
measured. Therefore, to fit these experimental data, the 
model depicted in Fig. [1] should be modified to include 
the bead detachment processes [see Fig. Elb)]. 

For the model depicted in Fig. EJb), the formulations 
of mean growth velocity V , mean growth and shortening 
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speeds Vg and Vg are the same as in Eqs. ([3]) and ^ 
^ . In the fohowing, we wiU get the expression of mean 
hfetime of the bead on MTs. Let Ti,T2,Ti, ,T2' be the 
mean first passage times (MFPTs) of a bead to detach- 
ment with initial sub-states_J_,__2^ 1' and 2' respectively, 
then Ti,T2,Tv,T2' satisfy 
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To 



1 



fci 



^1 ~l~ ~\~ kc 

1 



ki + ka + kc 
ks 

kr + k-i + k^ 



-To 



-Tv 



ki + ka + kc 



1 



1 



Ty 



(7) 



k2 ki 
Then the mean lifetime can be obtained as follows 

T = piTi + P2T2 + PlTy + P2T2' , (8) 



where pi,p2, pi, P2 can be obtained by formulations in 
Eq. @. 

In assembly state, let Tai and Ta2 be the MFPTs to 
detachment of the bead initially at sub-states 1 and 2 
respectively, then Tai,Ta2 satisfy 



1 ki 1 

Tal = -, —, h — — Ta2, Ta2 = 1 V Tal- 

ki + ka ki + ka k2 



One can easily show 

rj. ki+k2 



T 



a2 



(9) 



(10) 



k2ka k2ka 

Therefore, the MFPT to detachment of the bead in as- 
sembly state is 



Ta = PlTal +P2Ta2, 

where the steady state probabilities 



ko 



Pi 



P2 



ki 



(11) 



(12) 



ki + k2 ki + k2 

are obtained from Eq. The mean detachment rate 

during assembly can then be obtained by Ka = 1/Ta = 
l/ipiTai -\-p2Ta2), i.e., 



Ka 



(fci -I- k2)k2ka 

{ki + k2Y + kika' 



(13) 



Similarly, the mean detachment rate during disassembly 
can be obtained as follows 



Kd=l/Td^l/{piTdi+p2Td2) 
_ (fcs -I- k4)kikd 
~ {ks + ki)^ + kskd' 



(14) 



with steady state probabilities pi = fc4/(fc3 -I- ^4), p2 
ks/iks + ki). 



Let Tci and Tc2 be the MFPTs of MT to catastrophe 
from sub-states 1 and 2 respectively, then T^i and Tf,2 
satisfy [see Fig. ^h)] 

^ + T^^^^^ Tc2 = ^+Tci. (15) 

fci + kc ki + kc k2 

The mean rate of catastrophe can be obtained by Kc ~ 
1/Tc = l/(piTci +p2Tc2)- The explicit expression can be 
obtained by replace ka with kc in Eq. ()13|) 



Kc 



{ki + k2)k2kc 
(fci + k2f + fcifcc 

Similarly, the mean rate of rescue is 
(fca -|- k4)kikr 



Kr = 



(fca + fc4)2 + k^kr 



(16) 



(17) 



B. 



Force and GTP-tubulin concentration 
dependence of the transition rates 



From the experimental data in or see Fig. one 
sees some transition rates in our model should depend on 
the external force. Since the processes 1 — 2 and 2' <— 1' 
are accomplished by binding tubulin unit to and releasing 
tubulin unit from the tip of PF [see Fig. [l]and[2lb)], we 
assume that fci and fcs are force independent. Similar 
as the methods demonstrated in the models of molecular 
motors 26[ 38 1 and models for adhesive of cells to cells 
[3^, the external force F dependence of rates fc2, fc4, fca, 
kd, kr, kc are assumed to be the following forms 



fc; = fcOe^^*'/'=«'^, / = 2,4, a, d,r,c. 



(18) 



Hereafter, the external froce F is positive if it points to 
the direction of MT growth. 

Meanwhile, the rate fci should depend on the concen- 
tration of free GTP-tubulin in solution. Similar as the 
method in [26| , we simply assume fci — k^ [Tubulin] . 

C. Critical forces of MT growth 

For the sake of the better understanding of external 
force F dependent properties of MTs and the experimen- 
tal verification of the two-state model, in the following, 
we will define altogether eleven critical forces. Corre- 
sponding numerical results will be presented in the next 
section. 

(1) Critical Force Fci: under which Vg{Fci) = Vs(Fci), 
i.e. the average speeds of assembly and disassembly are 
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the same. From Eqs. ([5]) ^ one sees Fd satisfies 

kik2{Fcl)[h + k4{Fcl)] = k3k4{Fcl)[ki + k2{Fcl)]. 

(19) 

(2) Critical Force Fc2- under which the mean veloc- 
ity of MT growth is vanished. Formulation (jS]) gives 
kikr{Fc2) = ks,kc{Fc2), i.e. 



kBT 



ksT 



In 



In 



(20) 



{5r-5c)L fcOfcO [Tubulin]' 



(3) Critical Force Fc^: under which pi-|-p2 = Pi+P2, i-e., 
the probabilities that MTs in assembly and disassembly 
states are the same. From expressions in Eq. ([5]), one 
easily sees Fc^ satisfies 



kr{F,s)k4{F,^)[ki + k2{F,^)] 

^kc{F,3)k2{Fc3)[k3 + k4{F,3)]. 



(21) 



(4) Critical Force Fd : under which the detachment rates 
during assembly and disassembly states are the same. In 
view of formulations (fT3|) and ([T4|) . one can get Fca by 
Ka{F,i) ^ Kd{F,i). 

(5) Critical Force Fc^: under which the mean dwell times 
of MT in assembly and disassembly states are the same. 

Let Tgi and Tg2 be the MFPTs of bead to detachment 
or catastrophe of MT with initial sub-states 1 and 2 re- 
spectively, then Tgi,Tg2 satisfy (see Fig. [2](b) and Refs. 

silii) 



1 



fci 



fci 
k2 

Its solution is 
fci- 



-Tgl. 



ki 



-T, 



g2, 



(22) 



k2{ka 



T, 
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ki + ka + kc + k2 

k2{ka + kc) 



(23) 



The mean dwell time of MT in assembly (or growth) state 
is then 



Tg = piTgi -I- P2T, 



(fci 



ki{ka 



/C2(fcl + k2){ka + kc) 



(24) 



Similarly, the mean dwell time of MT in disassembly (or 
shortening) state can be obtained as follows 



{k-j + kjf + k-i{kd + kr) 

kiiks + ki){kd + K) 



(25) 



The critical force F^^ can then be obtained by Tg{Fc^) 
Ts{Fa). 



(6) Critical Force Fcq: under which the mean lifetime 
of the bead on MT attains its maximum, i.e. T[Fcq) — 
maxi? T{F) with T given by formulation ([8|). 

(7) Critical Force Fc^: under which the mean growth 
length of MT attains its maximum. The mean growth 
length of MT can be obtained by 1+ = VT with V, T 
satisfy formulations ([3]) and ([8]) respectively. 

(8) Critical Force Fcs- under which the mean shortening 
length of MT attains its maximum. The mean shortening 
length of MT can be obtained by I- = -VT with V,T 
satisfy formulations ([3]) and ([8]) respectively. 

(9) Critical Force Fcg: the rates of catastrophe and res- 
cue are the same, i.e. Kc{Fcg) — Kr{Fcg) [see Eqs. (fT6|) 
and (|17[) ]. Under critical force Fcq, the average switch 
time between growth and shortening, i.e. 1/Kc and 
1/Kr, are the same. It is to say that the mean dura- 
tion for each growth and each shortening period are the 
same. 

(10) Critical Force Fciq: under which V^Tg = VsTs- Here 
VgTg =: Zg is the mean growth length before bead detach- 
ment or catastrophe, and VsT^ —'. Is is the mean short- 
ening length before bead detachment or rescue. The for- 
mulations of Vg,Vs and Tg,Ts are in Eqs. (O ^ and 

(11) Critical Force Fdi- under which Vg/Kc — Vs/Kr- 
Here Vg/Kc =: Ig is the mean growth length before catas- 
trophe, and Vs/Kr — I* is the mean shortening length 
before rescue. The formulations of K^ Kr are in Eqs. 
CH) dlTl). 

It needs to be clarified that, the definitions for 
Fci, Fc2, Fc3, FcQ, Fell are unrelated to bead detachment, 
but the definitions for Fc4, Fc5, Fee, Fcj, Fcg, Fdo do. 
Therefore the values of Fd, Fc2, Fc3, Fcg, Fdi obtained in 
this theoretical study can be verified by various experi- 
mental methods as in 13 , 30 - 32 , , but the values of 



Fci, Fc5, Fc6, Fc7, Fc8, Fdo can only be verified by similar 
experimental method as in 27[ . For the sake of conve- 
nience, and based on the above definitions and numerical 
calculations in Sec. Ill (see Figs. [7]and[8|), basic prop- 
erties of the eleven critical forces Fd are listed in Tab. 
|T1 Meanwhile, the main symbols used in this study are 
listed in Tab. HD 



III. RESULTS 

In order to discuss the properties of MT growth and 
shortening, the model parameters, i.e. ki,k3 and /e°, 5i 
for i = 2,A,a,d,r,c should be firstly obtained. By fit- 
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TABLE I: Basic properties of the critical forces as defined in 
Sec. II.C, see also Figs. [7|and[8l 



TABLE IIL Model parameters obtained by fitting the expres- 
sions in (O ® ([HI) (HI) iMl) ([161) (in|) to the experimental data 



i 


F < F^^ 


F = F,, 


F > Fcr 


mainly measured in (27|]. In the fitting, 


ksT = 4.12 pN-m 


1 
2 


Vg < Vs 

V <Q 


Vg = Vs 

V ^0 


Vg > Vs 

V >Q 


and ettective step size L = U.bio nm are used yy, [4^ 
The fitting results are plotted in Fig. [S] 
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Pi + P2 < Pl + P2 


Pi + P2 = pi + P2 


Pi +P2 > Pl + P2 


Parameter 


value 


Parameter Value 


4 


Ka < Kd 


Ka = Kd 


Ka > Kd 


k\ 5. 


B X 10^ s'VM"^ 




9.3 s-i 


c: 



rp \ T-T 

lg> la 


Ig - -Ls 


Ig < Is 


kz 


1.0 X 10** s-i 


kl 


571.6 s-i 
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T < max p T 


T = maxF T 


T < maxF T 




7.4 X 10^ 




1.8 X 10^ s-i 


7 


1+ < maxF 1+ 


1+ = maxF 1+ 


l-f- < maxF 1+ 


l^a 


41.9 s~^ 




1.4 X 10'' s"^ 


8 


/_ < maxi? I- 


I- — maxi? I- 


< maxF l- 




0.68 


5s 


-2.88 


9 


Kr < Kc 


Kr = Kc 


Kr > Kc 




3.71 


5c 


-2.96 


10 


Ig < Is 


Ig Is 


^g ^ 


5a 


1.77 


5d 


0.33 


11 






''g ^ ^s 











TABLE IL The main symbols and their expressions (or defi- 
nitions) used in this study. 



Symbol 


Biophysical meaning 


Definitions 


V 


mean velocity of MTs 


Eq. ^ 


Vg 


growth speed of MTs 


Eq. ([U 


Vs 


shortening speed of MTs 


Eq. ([SI 


T 


mean lifetime of bead 


Eq. ^ 


Ka 


bead detachment rate (assembly) 


Eq. ^ 


Kd 


bead detachment rate (disassembly) 


Eq. ([HI) 


Ta 


time to detachment (assembly ) 


l/Ka 


Td 


time to detachment (assembly ) 


1/Kd 


Kc 


catastrophe rate 


Eq. ^ 


Kr 


rescue rate 


Eq. ([TT]) 


Pi, pi 


probability 


Eq. ^ 


Tg 


mean growth time 


Eq. ^ 


Ts 


mean shortening time 


Eq. ^ 


Ig 5 


Ig — VgFg ^ Is — VsTs 


see Fcio 


^g 1 


/; = Vg/Kc, i: = Vs/Kr 


see Fell 




1+ = VT, I- = -VT 


see Fc7,8 


ki 


rate constants [Fig. [2j2)] 


Eq. ([TH) 



ting the expressions of Vg,Vs,T, Ka, Kd, Kc, Kr, which 
are given in Eqs. (O (O ([U (dH) ^ ^ ^ respec- 
tively, to the experimental data mainly measured in [27| , 
these parameter values are obtained (see Fig. [3] and Tab. 
1111) the fitting methods are discussed in 46']). The data 
corresponding to zero external force in Figs[31[b) and[3Jc) 
[the two black dots on vertical axis] are obtained by fit- 
ting the corresponding measurement in [s^ with a con- 
stant [see the two lines in Figs. lUJa) andlUJb)], since as 
implied by our model, the rates of MT shortening and res- 
cue are independent of GTP-tubulin concentration. All 
the following calculations will be based on the parameters 
listed in Tab. IIIIl The curve in Fig. HJb) is the theoret- 
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FIG. 3: Theoretical results of the two-state model [see Fig. 
^h)] and experimental data obtained by Akiyoshi et al [27J: 
The detachment rates are obtained by formulations (|13|l (I14|) . 
the speeds are obtained by formulations ([Sj ([6}, the switch 
rates are obtained by formulations p6|) p7|) . and the lifetime 
is obtained by formulation ([8]). The model parameters used 
in the theoretical calculations are listed in Tab. [ml The 
two black dots on vertical axis of (b) and (c) are obtained by 
averaging the data in [30| [see the lines in Fig. |4ja) and Fig. 
Hb)]. 

ical prediction of GTP-tubulin concentration dependent 
catastrophe rate Kc by formulation (1161) . Compared with 
the experimental data measured in [sdls^l , these predic- 
tion looks satisfactory 47 [. 

From Fig. [Sfa), one can see that, the MT is mainly in 
assembly state. Further calculations indicate that the ra- 
tio of probabilities in assembly state to disassembly state, 
i.e. {pi +P2)/{pi + P2), increases exponentially with ex- 
ternal force F [see Fig. |6l[a)]. In experiments of Akiyoshi 
at al 27[, the external force F is applied to MT through 
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FIG. 4: GTP-tubulin concentration dependent data measured 
by Walker et al [s^]: (a) Shortening speed of MT and their 
average value, (b) Switch rates of MT between assembly and 
disassembly, where the curve is obtained by our theoretical 
model using the parameter listed in Tab. IIIII (see formulation 
(HSJ with fci = fe? [Tubulin]), the solid squares are experimen- 
tal data from 13211. 



a bead attached to its growing tip. Fig. ^h) indicates 
that, for F < 16 pN, the mean dwell time of MT in as- 
sembly state before bead detachment is larger than that 
in disassembly state. Although the MT is mainly in as- 
sembly state, its mean growth velocity is negative under 
small external force [Fig. [S^c)], since for such cases, the 
shortening speed is greatly larger than the growth speed 
[see Fig. [2Ib)]. But, Fig. [SJc) indicates the mean veloc- 
ity of MT growth always increases with external force. 
Similar as the mean growth velocity, the mean growth 
length of MT before bead detachment might be negative 
[i.e. MT shortens its length in long time average, see Fig. 
Eljd)], though the MT spends most of its time in assem- 
bly state [Fig. [SJb)]. Similar as the mean lifetime [Fig. 
[3]Jd)], the mean growth length of MT also has a global 
maximum for external force [Fig. [5][d)]. As we have 
mentioned in the Introduction, the chromosome segrega- 
tion is accomplished by the tensile force generated during 
MTs disassembly. Fig. [S] tells us the critical force of one 
MT disassembly is about 1.2 pN under the present exper- 
imental environment 



27|. In Fig. El^b), the mean growth 
length lg,lg and mean shortening length ls,ll which are 
given in the definitions of critical force Fdo, Fdi are also 
plotted as functions of external force. One can easily see 
that Ig < I*, and Ig < I* since the mean dwell time of 
MT in assembly state Tg < 1/Kc and mean dwell time 
in disassembly state T, < 1/Kr- But for large external 
force. Is ~ II since, for such cases, MTs leave disassembly 
state mainly by rescue. 

Since the assembly of MT depends on free GTP- 
tubulin concentration (in our model, the simple relation 
ki — /cj' [Tubulin] is used, and the disassembly process 
is assumed to be independent of GTP-tubulin concen- 
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FIG. 5: Properties of MT growth and shortening obtained 
by the two-state model [see Fig. |2{b)] with model parame- 
ters listed in Tab. IIIII (a) Under external force, the MT is 
mainly in assembly state, both the probabilities pi and p2 of 
MT in assembly sub-states 1 and 2 increase with force but 
with p2 S> pi. During disassembly state, the MT is mainly 
in sub-state 2', p2 > pi- Here pi,p2, pi, p2 are calculated by 
formulations in ([2]). (b) The dwell time Tg of MT in assembly 
state is always larger than that in disassembly state (denoted 
by Ts) for external force less than 16 pN [see formulations (|24|) 
(|25|) for Tg and Ts]. Similar as the mean lifetime of bead at- 
tachment to MT [see Fig. El^d)], both Tg and Ts increase firstly 
and then decrease with external force, (c) The mean velocity 
of MT growth [see formulation (|3}] increases with external 
force monotonically, where the negative velocity means MT 
shortens its length in long time average, though the curves in 
(a) and (b) imply that the MT is mainly in assembly state, 
(d) The mean growth length before bead detachment increases 
firstly and then deceases with external force. Here the mean 
growth length is obtained by mean growth velocity of MT 
multiplied by mean lifetime of the bead, i.e. VT [see formu- 
lation ^ for V and formulation (O for T]. 



tration, which can be verified by the data in (3o| . see 
Fig. HI, the eleven critical forces defined in the pre- 
vious section also depend on GTP-tubulin concentra- 
tion. For convenience, in our calculations (the results 
are plotted in Figs. [7]and|S]), [Tubulin] =1 means the 
free GTP-tubulin concentration is the same as the one 
used by Akiyoshi et al 0. From Figs. [7] and H 
one can see, the critical forces Fd for i ~ 4, 5, 6 in- 
crease, but others decrease with GTP-tubulin concen- 
tration [Tubulin]. For high GTP-tubulin concentration, 
Fc2 ~ ^cii and w Fcg since for kc <C fci, equations 
Vg/Kc — Vs/Kr and Kc = can be well approximated 
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FIG. 6: (a) The ratio of probability pi + p2 that MT in as- 
sembly state to probability pi + p2 that MT in disassembly 
state increases exponentially with the external force, and un- 
der positive external force, the MT mainly stays in assem- 
bly state, although the assembly speed might be much lower 
compared with the disassembly speed [see Fig. El^b)]. (b) The 
mean growth length Ig, I* and shortening length 1^, of MT in 
one assembly and disassembly period. The difference between 
Ig, Is and Ig, II is that, in the calculation of 1^,1^, the bead at- 
tached to the tip of MT, through which the external force 
is applied to MT, is assumed to keep attachment to MT, or 
the external force just applied by other methods [jol - li^ . |4]]| . 
so the MT can only leave assembly state by catastrophe and 
leave disassembly state by rescue. 




[Tubulin] 



[Tubulin] 



FIG. 7: Critical forces as defined in Sec. II. C, in which [Tubu- 
lin] =1 means the concentration of GTP-tubulin is the same 
as the one used in the experiments of Akiyoshi et al [i^j . For 
better understand the curves for Fd, see Tab. [D 



by kiKr = k^kc and krki{ki + ^2) ~ kck2{k^ + k/i) [4 
Since the force distribution factors 5g > Q, 5^ < Q (see 
Tab. imp , from Eqs. ([5]) ([6]) one can easily show that 
the growth speed Vg increases but the shortening speed 
Vs decreases with external force F. Therefore, Vg < Vg 
a F < Fci (see Tab. Eqs. (P ^ also indicate that 
the growth speed Vg increases with but the shortening 
speed Vs is independent of GTP-tubulin concentration 
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FIG. 8: Plots of cntical forces Fd for i = 1, 2, 3, 6, 9, 11 under 
low GTP-tubulin concentration. The meaning of [Tubulin] is 
the same as described in the caption of Fig. [T] 



[Tubulin]. Therefore, the critical force Fd decreases with 
GTP-tubulin concentration [Tubulin] [see Fig. [SJ^a)]. But 
for high [Tubulin] , critical force Fd is almost a constant 
[see Fig. [Tja)] since, for saturating concentration, the 
growth speed Vg tends to a constant [see Eq. ([5]) and 
Fig. El^a)]. The decrease of critical force Fc2 with [Tubu- 
lin] can be easily seen from expression (|20|) [see Fig. El^b)]. 
The decrease of critical forces Fd , Fc2 implies that low 
GTP-tubulin concentration might be helpful to chromo- 
some segregation. From expressions in ([2|) one can verify 

(Pl + V2)I{P\ + Pi) = krki^^x + fc2]/^cfc2[^3 + ^4]- So 

(pi +P2)/(pi + P2) increases linearly with [Tubulin] [see 
Fig. [TOT a)]. At the same time, (5^ + (5s > 0, (5g > and 
bg + ba < 0,(5s < (see Tab. imply (pi -f P2)/(pi +P2) 
also increases with external force F [see Fig. [^a)]. 
Therefore, the critical force F^z decreases with [Tubulin] 
[see Fig. [H^c)]. 

Since the detachment rate Ka increases and detach- 
ment rate decreases with external force F [see Fig. 
[3l[a)], and Ka increases with but K^^ is independent of 
[Tubulin] [see Eqs. (HH)], the critical force F^ in- 

creases with [TubuUn] [see Fig. [DJa)]. The increase of 
critical force Fc5 indicates MTs will spend more time in 
assembly state at high GTP-tubulin concentration [see 
Tab. Uand Figs. Elja) and [51(b)]. The increase of critical 
force Fee [see Fig. [8l[d)] implies, the peak of the lifetime- 
force curve as plotted in Fig. [3l[d) will move rightwards as 
the increase of [Tubulin] , but with a upper bound around 
4 pN [see Figs. Hid) and[9l[d)]. Similarly, the decrease of 
critical force Fc7 [see Fig. [71[d)] means, the peak of the 
mean growth length-force curve will move leftwards as 
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the increase of [Tubulin], and with lower bound around 
4.44 pN. Finally, critical forces F^s, F^g, F^o, F^i all 
decrease with [Tubulin]. It may need to say that, in Ref. 
27| . only experimental data for positive force cases are 
measured, and similar experimental methods as used in 
Refs. 32, |4l| might be employed to apply negative force 



to MTs. At the same time, the mechanism of MT growth 
and shortening under negative external force cases might 
be completely different from that under positive external 
force cases, so for the results of critical forces plotted in 
Fig. [7] which have negative values, experimental verifica- 
tion should be firstly done before further analysis. 

To better understand the GTP-tubulin concentration 
[Tubulin] dependent properties of MT assembly and dis- 
assembly, more figures are plotted in Figs. [9]and[T0l One 
can see that the mean lifetime T, ratio {pi+P2)/{pi+P2), 
and mean growth length lg,lg all increase linearly with 
[Tubulin] (from the corresponding formulations, one can 
easily see that the mean shortening speed Vg, and mean 
shortening length Is, I* are all independent of [Tubulin]). 
The mean velocity V and mean growth speed Vg also 
increase with [Tubulin], but tend to a external force F 
dependent constant [one can verify this limit constant is 
fcjL — k2Lexp{FSg/kBT). For such cases, the MT stays 
mainly in sub-state 2, i.e., p2 ~ 1, see Fig. [SJa)]. The 
mean growth length VT does not change monotonically 
with external force [see Figs. El^a) andHl^c)] but increases 
with [Tubulin] for high GTP-tubulin concentration cases. 



IV. CONCLUDING REMARKS 

In this study, a two-state mechanochemical model of 
microtubulin (MT) growth and shortening is presented. 
In assembly (growth) state, one GTP-tubulin will attach 
to the growing tip of the protofilament (PF) firstly and 
then, after the hydrolysis of GTP in the penultimate PF 
unit, the curved PF segment is slightly straightened with 
one new PF unit lying into the MT cylinder surface. In 
disassembly (shortening) state, one tubulin unit will de- 
tach from the tip of PF, and then the GDP (or GDP+Pi) 
capped tip segment of PF will be further curved with 
one new tubulin unit out of the MT surface (the phos- 
phate is assumed to be released simultaneously). The PF 
can switch between the assembly and disassembly states 
with external force dependent rates stochastically. Each 
assembly or disassembly process contributes to one step 
of growth or shortening of MT with step size L =0.615 




Force (pN) 



FIG. 9: GTP-tubulin concentration dependent properties of 
MT growth and shortening, (a) The mean growth veloc- 
ity V [see formulation Q] increases with GTP-tubulin con- 
centration. The plots also indicate V increases with exter- 
nal force [see Fig. [5{c)]. (b) and (d) The mean lifetime of 
bead attachment to MT increases with GTP-tubulin concen- 
tration, but increases first and then decreases with external 
force. For high GTP-tubulin concentration, the critical force 
Fed under which the mean lifetime gets its maximum is al- 
most a constant (about 4 pN), see also Fig. [TJd). (c) The 
mean growth length of MT before bead detachment does not 
change monotonically with external force and GTP-tubulin 
concentration, so there exists critical force under which the 
maximum is obtained. But for high GTP-tubulin concentra- 
tion, mean growth length increases with [Tubulin], (d) For 
any GTP-tubulin concentration, the mean lifetime does not 
change monotonically with external force. The optimal value 
of external force, under which the mean lifetime is maximum, 
increases with GTP-tubulin concentration, but is almost in- 
variable for large [Tubulin], (c) (d) Both the mean growth 
length and mean lifetime do not change monotonically with 
external force, so there exists optimal values under which the 
corresponding maximum is reached (see Fc6,Fc7 Fig. [7|. 



nm. This model can fit the recent experimental data 
measured by Akiyoshi et al [27[ well. 



From this model, interesting properties of MT growth 
and shortening are found: Under large external force or 
high GTP-tubulin concentration, the MT is mainly in as- 
sembly state; The mean lifetime of bead attachment to 
MT and mean growth length during this period (in ex- 
periments, the external force is applied to MT through a 
bead attached to the growing tip of MT) increase firstly 
and then decrease with the external force, but roughly 
speaking, they all increase with the GTP-tubulin concen- 
tration; The growth speed of MT increases with GTP- 
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FIG. 10: The ratio of probability pi + p2 that MT in as- 
sembly state to probability pi + p2 that MT in disassembly 
state, and the mean growth length lg,lg in each growth pe- 
riod all increase linearly with GTP-tubulin concentration (a) 
(c) (d) [see Tab. |TT]for definitions]. The growth speed Vg of 
MT increases with [Tubulin] but tends to a limit constant for 
saturated concentration [see Eq. ^ for the formulation of 

Vg]. 



tubulin concentration but has an external force depen- 
dent limit. For the sake of experimental verification, al- 
together eleven critical forces are defined, including the 
force under which the mean lifetime or mean growth 



length reach its maximum, the mean assembly speed is 
equal to the mean disassembly speed, the probabilities of 
MT in assembly and disassembly states are equal to each 
other, the detachment rates of bead during assembly and 
disassembly states are the same, the mean dwell times in 
assembly and disassembly states are the same, the mean 
growth velocity of MT is vanished, etc. Almost all of 
the above critical forces decrease with the GTP-tubulin 
concentration, since high GTP-tubulin concentration is 
favorable for MT growth and under low GTP-tubulin 
concentration, MT will shortens its length in average. 
Roughly speaking, GTP-tubulin and external force are 
helpful to MT assembly, but there exists optimal values 
external force for the mean lifetime of bead on MT and 
mean growth length of MT. 
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tion in [23] is different from that in [41] ]. In this study, 
the force direction definition is the same as in i.e., 
the force is positive if it points to the MT growth direc- 
tion), so the corresponding force distribution factor 5g 
[see Eq. (|18p ] should be positive since the growth speed 
Vg — kik2L/{ki + ^2) [see Eq. ([5])]. Consequently, the 
probability pi — fc2/(fci -1-^2) that MT in state 1 [see 
Eq. (I12|l ] increases but the probability p2 = fei/(fci -I- fe2) 
that MT in state 2 decreases with external force, i.e., as 
the increase of external force, the assembly MT would 
more like to stay in state 1. Meanwhile, from the exper- 
imental data in one sees the detachment rate from 
assembly state increases with external force. Therefore, 
the more reasonable choice is to assume that the bead 
can only detach from state 1 but not state 2. Through 
similar discussion, one also can see that it is more rea- 
sonable to assume that, in disassembly state, the bead 
can only detach from state 1'. At the same time, the 
experimental data in 30] imply the catastrophe rate 
decreases with GTP-tubulin concentration [Tubulin] [or 
see Fig. HJb)]. Since ki = fc? [Tubulin] , the probability 
pi = ^2/(^1 + fe) decreases with [Tubulin], but the prob- 
ability p2 = k\/{k\ -\- k^) increases with [Tubulin]. This 
is why we assume the catastrophe takes place at state 1. 

[46] In our fitting, we firstly get the parameters fcj, k2,Sg and 
kg, k1,Ss by fitting formulations ((5)1 and ([Sjl to the experi- 
mental data of growth and shortening speeds respectively 
[see Fig. |3jb)], and then get k^,Sr and kc,3c by fitting 
formulations (|16|l and (|17|) to the catastrophe and res- 
cue rates [see Fig. |3][c)], k^l,Sa and k'^jSd are determined 
by fitting formulations (I13p and (I14p to the correspond- 
ing data plotted in Fig. [Sfa). Finally all the parameters 
are slightly adjusted according to the experimental data 
about the mean lifetime of bead attachment to MT [see 
formulation ^ and Fig. [3jd)] . All the fitting are done by 
the nonlinear least square program IsqnonUn in Matlab. 
In each fitting. We randomly choose 1000 initial values 
of the parameters and adopt the parameter values which 
fit the experimental data best. 

[47] The parameter values listed in Tab. Illll do not fit well to 
the GTP-tubulin concentration dependent growth speed 
of MTs obtained in [sol . |33 ] , since the corresponding data 
in [iol . [3^ are much different from that in [23] . Without 
external force, but under similar GTP and tubulin con- 
centration, the growth speed of MT measured in [i^l is 
about 43 nm/s, and about 20 nm/s in [s^], but it is only 
about 5 nm/s in [23]. In this study, we get the param- 
eter values mainly based on the data measured in EtI . 
One reason is that, from our model, if the GTP-tubulin 
concentration is nonzero, the growth speed of MT will 
always positive [see formulation (O, Vg > if fci >_0l. 
However, this might not be true for the data in [jol. |32|. 
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So, it might be impossible to get a believable fitting pa- 
rameters for formulation ([5]) from data in smce 
the data in [jol . Is^ ] cannot be described by a formulation 
like ([5]). In [4l|], the velocity- force data are measured un- 
der tubulin concentration 25/iM. However, the zero force 
growth speed obtained there is about 20 nm/s, which 
is also much larger than that obtained in [23] • Conse- 
quently, the theoretical results based on the parameter 
values listed in Tab. IIIII do not fit well to their data ei- 
ther. One can verify that the velocity-force data in [Zlj ] 
can be well described by formulation ((5)l but with param- 
eters fc? = 2.99 s-VM"\ fc§ = 53.17 s"^ and 5g = 4.13. 
The difference among these experimental data might due 
to the differences of experimental techniques, methods or 
materials. 

[48] If the simple model depicted in Fig. [2ja) is employed to 



describe the dynamic properties of MT, then Fc2 = -Fcii 
and Fc3 = Fcq. The reason is as follows. At steady state, 
the probabilities p, p that MT in assembly and disassem- 
bly states are p — fc,./(fcc -I- fcr) and p = kc/(kc + kr) 
respectively. So the mean growth velocity of MT is V = 
k\p — k2p = (fcifcr — k2kc)/(kc H- kr). Then the criti- 
cal force Fc2 satisfies ki{Fc2)kr{Fc2) — k2{Fc2)kc{Fc2)- 
Meanwhile, I* — Vg/Kc ~ Vg/kc = k\L/kc and I* = 
Vs/Kr = Vs/kr = k2L/kr, SO l*g{Fcii) = l*s(Fcii) is 
equivalent to fci (Fcii)fcr(-Fcii) = k2(Fcii)kc{Fcii) which 
means Fc2 ~ Fdi. At the same time, p — p is equivalent 
to kc — kr, so Fc3 — Fc9. But for our model as depicted 
in Fig. Mh), Fc2 / -Fell and F^-j. ^ F^g [see Figs. Mh) 
and [He)]. 



